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Abstract 

Self-dual vortex solutions are studied in detail in the generalized abelian Higgs 
model with independent Chern-Simons interaction. For special choices of couplings, 
it reduces to a Maxwell-Higgs model with two scalar fields, a Chern-Simons-Higgs 
model with two scalar fields, or other new models. We investigate the properties of 
the static solutions and perform detailed numerical analyses. For the Chern-Simons- 
Higgs model with two scalar fields in an asymmetric phase, we prove the existence 
of multisoliton solutions which can be viewed as hybrids of Chern-Simons vortices 
and CP^ lumps. We also discuss solutions in a symmetric phase with the help of the 
corresponding exact solutions in its nonrelativistic limit. The model interpolating all 
three models — Maxwell-Higgs, Chern-Simons-Higgs, and CP^ models — is discussed 
briefly. Finally we study the possibility of vortex solutions with half-integer vorticity 
in the special case of the model. Numerical results are negative. 



I. INTRODUCTION 



In the last few years there has been much interest in the (2+l)-dimensional 
abehan Higgs models and the CP^ model with a Chern-Simons term in the gauge 
field action, partly because of their possible relevance in condensed matter physics 
[H]. These models allow classical vortexlike solutions and the analysis of them has 
become quite active recently 0, ^. Here particularly interesting are the so-called 
self-dual systems. Since the discovery of the self-dual charged vortex solutions in rel- 
ativistic pure Chern-Simons-Higgs system [Q, several other self-dual systems includ- 
ing the Chern-Simons term are now known |^. Of some interest among these 
is the generalized self-dual abelian Higgs model with independent Chern-Simons in- 
teraction, discussed in Refs. 0. It has two abelian gauge fields, one of which 
has only the Maxwell term as its kinetic term and the other only the Chern-Simons 
term. It may be considered as a realistic model if one identifies the Maxwell field as 
the electromagnetic field and the Chern-Simons field as an effective field arising in 
a condensed matter system. Besides, as we will see below, it contains a linearized 
version of the CP^ model with or without a Chern-Simons term as special cases 
and therefore interpolates abelian Higgs systems and the CP^ model. The present 
paper is devoted to the detailed study of static solutions in this model, which was 
not given in Ref. 0. 

This model consists of two complex scalar fields, a neutral scalar field, and two 
gauge fields. Each complex scalar field couples to both gauge fields and so there 
are four gauge couplings in general. But, with some special choices made for these 
couplings, the model reduces to simpler systems. For example, for some particular 
choices it becomes a decoupled sum of self-dual Maxwell Higgs and self-dual 
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Chern-Simons Higgs systems. More interestingly, it is possible to obtain a system 
in which two scalar fields interact with a Maxwell or Chern-Simons gauge field (but 
not with both), while possessing an additional global SU{2) symmetry. This kind of 
model with a Maxwell gauge field has attracted much attention recently ||^ because 
it admits topological vortex solutions in spite of having a simply connected vacuum 
manifold, and a class of solutions obtained here interpolates ordinary Ginzburg- 
Landau vortices and CP^ lumps [§]. We consider the corresponding model with a 
Chern-Simons gauge field, and find that it has nontopological vortex solutions in a 
symmetric phase as well as topological vortices in an asymmetric phase (as in the 
Maxwell case). Again a class of solutions obtained in this case approximate CP^ 
lumps, and the existence of multivortex solutions can be shown by adapting the 
method of Wang |^. There are also more general self-dual systems based on the 
same set of fields but with no additional global symmetry. 

If both gauge fields couple to scalar fields nontrivially, a new kind of models 
are obtained 0. If one of the two complex scalar fields decouples and the other 
couples to both gauge fields, only a broken phase is possible and topological vortices 
exist. We find that the nature of the vortex solutions is very similar to those of the 
Ginzburg-Landau system, the Chern-Simons field playing a minor role only. We will 
briefiy discuss its extension to the two complex-scalar-field case in such a way that 
there is an additional global SU{2) symmetry. Complicated but interesting is the 
case (f) of Ref. [Q, where all fields couple to one another nontrivially. In ref 
it was noted that there might be solutions with a half-integer vorticity (as regards 
each of the two gauge fields). We will report the result of our numerical work, which 
is negative. 

This paper is organized as follows. In Sec. II, we briefiy review the model to be 
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considered by us, following Ref. 0. Section III is the major part of the paper. We 
study static solutions for various interesting cases of this model, based on analytical 
and numerical means. We summarize our findings in Sec. IV. In the appendix, 
a f/(l) self-dual Maxwell-Chern-Simons Higgs system with additional global SU{2) 
symmetry is given. This model might be interesting because the solutions interpolate 
all three cases — Ginzburg-Landau, Chern-Simons, and CP^ solitons. 

II. THE MODEL 

We consider the (2+l)-dimensional system described by the following Lagrangian 

i 

C = —F^'F,, + l/te^'^^S.a, + \D,x? + I^^m^P + \{d^Nf - U{x,i^,N), (1) 
where D^x = [df, - «(eia^ + e2A^)]x, D^ip = [d^ - iie^a^ + eiA^)]ip, and 
Uix,i^,N) = he2\x\' + cTie,\^\^-uY + ^iei\x\' + a,eM^-a2—N-vy 

+ %mei\x? + aieslxP - ^2^N - v'f. (2) 

Here, ai and a2 may assume the values ±1 independently, x and ip denote complex 
scalar fields with couplings to both Maxwell and Chern-Simons gauge fields (i.e., 
and a^), while is a neutral scalar. The specific form of the potential (0) leads to 
a self-dual system for general (T2 = ±1, and only the case of di = cr2 = +1 was 
considered in Ref. P|. We shall now briefly summarize the results of Ref. and 
extend them to generally allowed values of ai, (T2. The variation of Oq and Aq leads 
to two Gauss laws 

«:/i2 = eiJj + esJj, (3) 
= e2Jj + e44, (4) 
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where /^,, = d^ay-d^a^, = -i{x*D'^X-Df^x*x), and = -i{tp* Df^tp - D^'4)*tp) . 
For finite energy configuration, Eqs. @ and @ imply relations among the magnetic 



flux = jd?xfi2 and charges Q^^ = jdPxJ^ and Q,/; = jd?xJ^, viz., 



= e2Qx + e4Qv 
There is no restriction on (^a = ^d?xFi2. The energy functional is 
E = l£x 



^ ^2 , 1 772 



(5) 
(6) 



+1 AxP + I A^P + ^(5oiV)2 + l(9.iV)2 + U 



After employing the Gauss laws and integrating by parts, one has the bound 



(7) 



(8) 



where the equality holds if and only if the following self-duality equations are satis- 
fied: 



{Di ± iD2)x = {Di ± mi A)^ = 0, 
i"i2±(e2|x|' + fXie4|^|'-M') =0, 
i^o ± <y2d;N = 0, 

Ax ± ^-x{e^\x? + ^ie,W - o^—N - v') = 0, 



A^ ± ^-^(eilxl' + t^iesl^P - a2— iV - v^) 

K 63 



0. 



(9) 



Note that if ai = — 1, the first equation shows that one of the two complex scalar 
fields is self-dual while the other is antiself-dual. Also, depending on the value of (T2, 
the sign in front of gets changed. However, each choice of (ai, 02) corresponds to 



a different potential as seen in Eq. and all signs in Eq. are correlated in a 
given theory. Each self-dual system for given ai and a2 is in fact not much different 
from the others. From Eqs. (|I|), (|^) and (|^), we see that the system with (Xi = — 1 
is transformed into that with ai = +1 under 63 — > —63, 64 — —64 and ip ip*, and 
into the system with a2 = —1 hj N ^ —N. Therefore, we put (Xi = cr2 = +1 from 
now on, except at the end of Sec. IIIA where we will give some further comments. 

We are interested in static solitonhke solutions of Eq. (|). Since there are four 
different gauge couplings, the model has a rich variety and the nature of solutions 
depends on the choice of these couplings. It should be thus desirable to consider 
separately various cases corresponding to special choices of couplings. For certain 
choices of couplings some of known self-dual systems, including the model of some 
recent interest [§] , are recovered. For other choices we encounter models which have 
not been studied in detail so far. 

III. SPECIAL CASES 

A. 61 = 63 = 

If the couplings associated with a^, i.e., ei and 63, are identically zero, then 
decouples and we have a self-dual Maxwell-Higgs system with two complex scalar 
fields and a neutral scalar. This is a generalization of the well-known self-dual 
Ginzburg-Landau model [0. Furthermore, from Eq. (^), the equation for be- 
comes V^A^ = 2(e2|xP + e4|^/'p)A^, which means that N may be taken to be zero 
identically. Here the potential is effectively 

U = lie2\x\' + eM'-^^Y■ (10) 
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If one further restricts oneself to the case 62 = 64, this model will possess additional 
global SU{2) symmetry and become in fact identical to the one recently considered 
by Vachaspati et al, Hindmarsh, and Gibbons et al. [If. For general 62 and 64, 
there is no additional global symmetry and we study this general case below. 
From Eqs. (0)and (p!OD , finiteness of energy requires that 

r ^ 00 : e2|xP + e4|V^P — ^ (11) 

as well as 

r^oo: Ax — 0, — 0. (12) 

Eq. (|ll|) tell us that on the circle at infinity, scalar fields must lie on an ellipsoid[] 
in the 4-dimensional internal space; then, Eq. (p!^ ) further restricts the scalar fields 
to be at most a pure phase there. Hence even if the vacuum manifold is simply 
connected, there can be vortex solutions characterized by a winding number 00. 

Now we consider the self-duality equations (P). On points where are nonzero, 
the first equation may be rewritten as 

i{Ai+iA2) = —{di+id2)\nx = — {di + id2) In ip, 

62 64 

where we chose the upper sign. Then 

{d, + td2)ln^ = 0, 

X 

so that 

^/,e2/e4 

w{z) ^ ^ (13) 

X 

^For the moment, we consider the case 62 > 0, 64 > only. 

^In fact this phenomenon has been known since seventies even if such model has not been 
explicitly constructed pO[ |. 
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is locally analytic in z = x + iy. Denoting Zr{r = 1, . . . , rii) and z'j.{r = 1, . . . , to 



be the zeros of x and ^ respectively, we may then write [g, ^ 



llr=ll^ ~ 



where h{z) has no zero or pole and is regular everywhere. Since x ^^id ip should 
approach vacuum values as r — > oo, h{z) is at most a constant, say g, by Liouville's 
theorem in complex variable theory. Note that w{z) is multivalued in general, but 
it is not problematic as long as x ^^id if) are single-valued functions. Now let 

where xo and ipo are some fixed constants satisfying e2|xoP + e4|V'oP = u^- If 62 = 64, 
all vacuums are equivalent because of global SU{2) symmetry and hence it suffices to 
choose any convenient values for xo and ipo, for example, V'o = and |xo| = ^fe^. 
However, if 62 7^ 64, various choices for the vacuum values are inequivalent and one 
must consider in particular the case with XoV'o 7^ 0. In the latter case, from Eq. 
(p!3D, we have wi^z) — > finite as r 00, which implies in view of Eq. (0) that 
62/64 = ni/n2. In other words, regular vortex solutions may exist only if 62/64 is 
rational and the ratio of the number of zeros for x and if), ni/n2, must have the 
same fixed value equal to 62/64. The fiux is then given by 

f , . 27mi 27rn2 , s 

$ = i dx'Ai = = -. (16) 

Jr=oo 62 64 

On the other hand, if one of the vacuum values, say ipQ, is zero we can have a 
totally different behaviors for the 62 7^ 64 case, which is absent for the case 62 = 64. 
From Eq. ([T3|), w —>■ as r 00 and this implies only the condition 62^1-2 < 6471,1. 



There is no need for 62/64 to be rational here and, except the above inequality, no 



condition on ni/n2, either. This can be understood once if one notices that the 
vacuum manifold for ip is trivial in this case. Now suppose that 

ih — > — as r — > oo. 

Then the flux will be, from Eq. (P), 

^ ^ 27mi ^ 27r(n2 + a) 

62 64 



or a is given by 



a = — ni-r22>0. (18) 

62 



There is a further non-trivial equation to be satisfied by |x| { oihy Elimi- 
nating gauge fields from the self-duality equations (^, we have in fact 

ni 



In - 262(621x1' + 64|V^r - = 4vr 5: 5(x - x,), (19) 

r=l 

where = {\w\\x\Y^^^'^ ■ If 62 = 64 = e, we can simplify this equation as 

V^r/ + 26^2(1 - e'') = In n 1-2 - -2-1' + 1^1' 11 k " . (20) 



\r=l r=l 

|2 I L;,|2-\ /„,2l 



where 77 = ln[e(|x| + IV^I )/"" ]■ Analytic solutions to Eq. ([191) or (pOf) are not 

available and we should resort to numerical analysis to go further. The simplest case 

of 62 = 64, with additional global SU{2) symmetry, gives rise to some interesting 

consequences, but they were extensively discussed already in Ref. 

For generic case, some numerical solutions have been studied by us for specific 

couplings and vorticities. In Fig. 1(a) we have plotted rotationally symmetric 

solutions with unit vorticity when the ratio p = — is equal to 1, or 2, while choosing 

the point |xoP = 12,62 on the ellipse of possible vacuum values. In addition, we 
■^However, see the next section. 



studied solutions for other n and p values and we found that the shapes are more 
or less the same for different p-values except that the vortex sizes become a little 
smaller for larger p. In Fig 1(b) we have plotted e^^^ for the case ni = 1 and n2 = 
with 62 = 64 = e. Note that the approximation e''/^ = 1 is excellent for |g| = 10 
and good even for |g| = 5. Therefore the solutions well approximate CP^ lumps for 
large as noted in Ref. For the number of free parameters in the general 

solution (which need not be rotationally symmetric), one can use the index theorem 
and this was already done in Ref. for the case with Xo'^Po 7^ 0. The number is 
equal to 2{ni + 712). 

Closing the discussions, we comment on the case with 6264 < 0. From Eqs. (p!3D 
and (0), we now find that 

1 ni n2 

- = = l[(z- ^,) Y[(z- 4)1^^/^^1/10, 

W r=l r=l 

where Hq is a constant. Then ni and ^2 should be zero because 1/w approaches 
a finite value as r — » 00. This implies that the energy E = u'^\^a\ is zero and 
there is no nontrivial solution satisfying the self-duality equations. However, as we 
mentioned in Sec. II, there is another self-dual system which corresponds to ai = —1 
in Eq. (H), and this system is essentially identical to that with 6264 > and Ci = +1. 

B. 62 = 64 = 

With 62 and 64 set to zero, the Maxwell U{1) field decouples from matters fields, 
and so does the neutral field A^. It is evidently the Chern-Simons counterpart of 
the model considered in the case A, and similar analysis can be used to find soliton 



* Readers may consult discussions given in the case B for the relevance of the CP^ model here. 
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solutions allowed for this case. Here we have the potential 

U=\{e\\x? + eM'){e,\x? + e,\^\'-v')\ (21) 

and, if ei = 63, there is an additional global SU{2) symmetry in the system. From 
the self-duality equations (^, we may conclude, just as in the case A (see Eq. (0)), 
that 

where q is an arbitrary complex constant and Zr{r = 1, . . . , rii) and zl.{r = 1, . . . , 77,2) 
are zeros of x and ip, respectively. In the present case, however, the situation is more 
complicated because there is also a symmetric phase as well as an asymmetric or 
broken phase; in the symmetric phase, we have an additional possibility of nontopo- 
logical soliton solutions. 

First, we concentrate on the topological soliton solutions allowed in the broken 
phase, for which the analysis is almost an exact parallel to that of the case A. 
For example, the statements below Eq. ([T^) apply also to this case with 62 and 
64 replaced by ei and 63, respectively. In particular, with ei = 63 = e, additional 
global SU (2) symmetry allows us to assume, without loss of generality, that 

(j)-(o)- (i-i^-g- (-) 

and we can follow Gibbons et al. @| closely. An immediate observation is that the 
function w{z) =ip /x should vanish at infinity and Eq. ( ^2]) reduces to 

where Pn{z) = l\r=i{^ ~~ ^r), and Qn{z) is a polynomial of z of order not larger than 
n — 1. Further, eliminating gauge field in Eq. (^ yields the following equation for 
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Ixl (note that 



4e 



In \^\'^ - ^{\^\^ + + _ _ \ 4^ ^ _ (25) 



r=l 



Choosing to work with dimensionless quantities, we make the replacements x 
fX/A/e, iIj — > vijj / y/e, 2v'^a^/K, and then Eq. (|25| ) reads 



In IxP - (IxP + m\{\x? + I^P) - 1] = 47r ^ 5(x - x,). (26) 

r=l 

To facihtate the analysis of Eq. (126|) , we introduce the quantity 

r7 = ln(|xr + |^r) (27) 



so that we can recast Eq. (E6]) as 



V'r/ - e''(e'' - 1) = p, 



(28) 



where we defined 

p = VMn(|P„p + |Q„|2)_ (29) 

[Note that In \z — Zrl"^ = 47r5(x — x^).] Since |xP + IV'P approaches 1 as r ^ cxd, 
we here expect that 77 — as r oo. Following ref. H], we define 



u 



1 = ln(|P„p + |Q„|2) _ ^ 1^,(1^ - + /i), 



(30) 



r=l 



where /i(> 0) is for the moment an arbitrary constant. Setting rj = ui + rj then gives 



V^fj-h + e"ie''(l - e^^e'') = 



(31) 



with h = A J2r=i p(l-2 ~ -2rP + p) ^- Eq- ( PTD corresponds to the variational problem 
associated with the functional 



A{fi)= j(fx ^|Vr7|2 + (e"i+*'-l)2 + 2/ir/ 



(32) 
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where fj should vanish sufficiently fast as r — ^ oo. With Qn = 0, this is in fact 
the action used by Wang ^ to prove the existence of multivortex solutions in the 
self-dual pure Chern-Simons model with one complex scalar field. His method can 
be adapted to the present case. Wang's proof for Qn = rests essentially on showing 
that A{fj) is a coercive functional and hence ^(^7) has a unique minimum. Other 
parts differing from the Qn = case having been checked already by Gibbons et al. 
[|] for sufficiently large /x, the only crucial step needed in establishing coercivity is 
to estimate the quantity jd?xe^'^^{e^ — 1)^. Wang proved for = that 

where C is a positive constant. However, it is not difficult to see that the same 
argument can be used also for the case 7^ and hence Eq. (|33[) holds even if 
7^ 0. We conclude that for every choice of polynomial P„ and Qn, a unique 
solution exists and the scalar fields may be reconstructed from r] via 



n 



V ^ / ^1 + V ^ / r=l k - Zr\ 

= I f ^" ^ e"^^ (34) 

^\Pn\'+\Qn\' V J 

up to gauge transformations. The gauge field Oj is then readily obtained using 
Eq. d^). The moduli space of solutions is just C^", the 4n-dimensional space 
parametrized by the complex coefficients specifying the polynomials P„ and Qn[|- 
We note that the equation of the form (|5^) has been given for the case A in Ref . ^ . 
The index theorem analysis of Lee, Min and Rim on the number of zero modes 

is also consistent with what we have shown rigorously here. 

^Explicitly, the free parameters can be identified with pk, qk [k = 0, 1, . . . , n — 1) when one 
writes Pn{z) = z" +p„-i2;""^ H hpiz +po and Qn{z) = (7„_iz""^ H h (?i + go- 
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A characteristic feature of Chern-Simons vortices is that, unhke the Ginzburg- 
Landau vortices, they carry nonzero angular momentum. The angular momentum 
is 



J = Jd^xe'^x\DoX*D,x + Do^*Do^ + c.c] 



cfxe'^x'- 



fl2 



x\^{aj-djcLTgx) + \'ip\'^{aj-djcLTgip) , (35) 



IxP + 1^ 

where we used the self-duality equations on the second expression. Now suppose we 
restrict ourselves to solutions with Pn{z) = z^^ , Qn{z) = qz"^^, {n2 <ni), where q is 
a constant, in Eq. ([3^). This corresponds to working with the ansatz 



(36) 



X = 






= 






tti = 


y,2 y 


ni+n2\ 
2 J 



where |^| and g are functions of r and |^|/|xl = 9^^^ • Then Eq. (^) becomes 
J = ^[/(oo)-g^(0)l + ^(n,-„,)jf *9'^snr^- '''' 

Obviously, g{0) = [rii + /2 from Eq. (0), while the finite energy condition 
requires that g{oo) = — |(^i — ^2)- Therefore, the first term on the right hand 
side of Eq. (|37|) is equal to —7cnnin2/e'^. But the second term cannot be explicitly 
evaluated unless ip vanishes identically. For the special cases corresponding toip = 0, 
we have 

J=— ^, (38) 



e 

which is in agreement with the result of Ref . @] . 
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Let us consider the case with ni = 1 and ^2 = in more detail. For this case we 
may write 

e^/^ (39) 



\ip J ^^2 |g|2 \^ q 

where we have set z — zi = re*^. It describes a vortex-hke structure centered at 
z = Zi, with size and orientation determined by the complex parameter q. Using 
Eq. p8|), we can determine the asymptotic behavior of 77: 

This shows that if g 7^ 0, scalar fields approach the vacuum values like 0(l/r^), 
contrary to the g = case for which the limit are approached exponentially. For 
g 7^ the magnetic field shows also a power fall-off at large r, being given by 

/.^^^i^M,o(l). (41) 



On the other hand, the behaviors for smaller r are found as 
= co + |(^^-co(co-l)^r2 + 0(r^), 

/12 = -^co(co-l)-^(l-2co)(^^-co(co-l))r2 + 0(/), (42) 

where cq is a constant (which is equal to 1 if g = 0). Note that the magnetic field 
is nonzero at the origin for g 7^ 0. Hence its behaviors are different from those of 
Ref. [Q at both short and large distances. Another interesting feature is that the 
solution approximates a CP^ lump when |g| ^ 1. This was noted in Ref. p| for the 
case A, and we can apply the same argument here. Namely, for |g| ^ 1, we have 
p = 4|gp/(r^ + |gp)^ ~ and rj ^ so that the scalar fields lie effectively on the 
vacuum manifold S^, 

(43) 



I ~ Jr2 + |g|2 V ^ 
14 



This represents a CP^ lump and we may say that the topological solitons of this 
model interpolate between Chern-Simons vortices and CP^ lumps. We have verified 
that the angular momentum J, given by the expression (|37|), approaches zero as we 
let g — s> oo. 

Some numerical results are shown in Fig. 2. In Fig 2(a) we plot rotation- 
ally symmetric solutions for the coupling ratios p = ^ = 1 and 2, while choosing 
|x(oo)p = IxoP = f^/2ei. As in the case A, sizes become smaller for larger p and the 
magnetic field gets confined in a narrow ring when p is large. The values of the mag- 
netic field are zero at the origin, which is a characteristic feature of Chern-Simons 
vortices 0|. In Fig. 2(b) we plot the solutions of the form (^) for |g| = 0.2,1,5, 
and 10. It is clear that the solutions approximate CP^ lumps when |q'|>5 just as in 
the case A. The magnetic field here shows interesting behaviors; for small |g| it is 
ring-shaped (as in Fig. 2(a)) even if it does not vanish at the origin, while it looks 
similar to that in Fig. 1(a) for large |g|. In this way, we see clearly that our solution 
indeed interpolates Chern-Simons vortices and CP^ lumps. 

Now we turn to nontopological soliton solutions allowed in symmetric phase. 
While Eq. ( ^2]) is still valid, there is no restriction on rii and n2 because the vacuum 
manifold is trivial. Let the asymptotic behaviors of scalar fields be such that 

, , 1 

101 - ^, (44) 

where a, P are positive constants. Then, by Eq. (^2]) , a and /? are not independent 
but related by 

63(^1 + a) = ei(n2 + /3). (45) 
This relation is due to the fact that there is only one gauge field in this case. The 
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flux is given by 

^ ^ 27r(ni + a) ^ 27r(n2 + (3) 
ei 63 

At present there is no rigorous existence proof of nontopological solitons in Cliern- 
Simons theory even with one scalar field, although their existence is strongly sup- 
ported by numerical analysis. The situation is similar for the present case. Here 
we expect that the number of zero modes be larger than 4n, as we have seen in the 
model with one scalar field |Q. 

We have done some analysis on nontopological solitons in the case of ci = 63 = e, 
assuming the form (|39|). Even if we do not know the exact solution, we can find 
approximate solutions for |g| ^ 1 by making use of the solutions for the correspond- 
ing non-relativistic self-dual model. Formally, in the non-relativistic limit, Eq. (|25|) 
reduces to 

n 

In \x\' + \x\' + I^P = 4vr 5: 5(x - x,), (47) 

r=l 

Other than the Liouville-type solutions obtained by setting |x| = {ipl, the present 



author found that Eq. (|47| ) admits a family of exact solutions ||12[ of the form 

/ X \ Vl2iPd,Q - Qd^P) ( P \ 



(48) 



\^ ) (|pp + |gp)3/2 \Q )^ 

where P, Q are polynomials of z sharing no common zeros and x^'s (see the right 
hand side of Eq. (|47|)) can be identified with zeros of P{PdzQ — Qd^P)- For our 
purpose, we take P{z) = e^^/^q^ z and Q{z) = e'^'^'^q^ so that 

f M = ( ] (49) 

J (r2 + |g|2)3/2 \ q J- ^ ' 

Eq. (|49|) will be a good approximate solution of the relativistic system when |g| ^ 1. 
We have found the asymptotic behaviors \x\ ~ r^^ and l^ipl ~ r^^, and so we identify 
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a = 2 and /5 = 3 for this solution. The gauge field Oj is given by Eq. ( p6D and g{r) 
has limiting values g{oo) = —5/2 and g{0) = —3/2. With these and using Eq. (|37D , 
we can calculate the angular momentum for this solution, 



J 



AlTK 


3ttk 








|g| > 1 


— 9 ' 



Therefore, in symmetric phase, we find a non-zero angular momentum even if |g| ^ 
1, contrary to the case in broken phase. Fig. 3(a) shows the plot of rotationally 
symmetric nontopological soliton solutions for a = 2.1 and 3 when p = 2. In Fig. 
3(b) we plot solutions of the form (p9|), for = 5 and 10. Note that for |g| = 10, it 
is well approximated by the exact solutions of nonrelativistic self-duality equations. 
It gives a = 2.056, which is not far from the nonrelativistic value a = 2. 

Closing this subsection we remark two things. First, this model may be viewed 
as a linearized version of the CP^ model with the Chern-Simons term. The po- 
tential has been chosen by demanding self-duality, and then the symmetric phase 
becomes degenerate with the asymmetric phase. In this sense, this model inter- 
polates two popular field-theoretic models useful for high-Tc superconductivity |1|, 
namely, the CP^ model and the Chern-Simons-Higgs model. Second, while the case 
A interpolates Ginzburg-Landau vortices and CP^ lumps, the case we just discussed 
interpolates Chern-Simons vortices and CP^ lumps (in broken phase). Then it is 
natural to expect that U{1) Maxwell- Chern-Simons theory with two complex 
scalar fields interpolates all three kinds of solitons. This model is briefiy discussed 
in the appendix. 
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C. 63 = 64 = 



In this case, the field ip decouples from the rest and the potential reduces to 

TT ^ / I |2 2\2 I ^1 I |2 / I |2 ^1 2 ^^^2 ,^\^ ,r--,\ 

U = ^[e2\x\ -u) +^\x\ ei X u , (51) 



where we have adjusted the field by a suitable constant shift. Assuming > 0, 
the vacuum configuration then clearly corresponds to 

= uVe2, AT = 0. (52) 

So only the broken phase is possible. This case is different from the previous two 
cases in that both Maxwell and Chern-Simons gauge fields couple to a single complex 
scalar x simultaneously. From Eqs. (|^) and (H), we note that 

K^a = exQx = 0- (53) 

Therefore, there is no charged vortex as noted in Ref. P|. If 62 = 0, this model 
reduces to the pure Chern-Simons model which is known to have charged vortices. 
Hence, 62 = is a singular point, i.e., the theory with 62 = has distinct behaviors 
from that with 62 7^ 0. 

We find it convenient to work with dimensionless quantities by making the re- 
placements 

^ ^=^, X ^ ^X, N ^ ^AT, 
V2e2-u Je^ kc^ 



72^ , f2 



All > \ I uA^. 



ei V ^2 

Then, after eliminating gauge fields, self-duality equations yield the following equa- 
tion for Ixl and A^ (here ^ = elu^ / K^e\): 



V^\Ti\x?-2i\x\Wx?-N-l)-\x? + l = 47r^5(x-x,), 

r=l 

V^N-\x\Wx?-N-l) = 0, (54) 



For these coupled equations, no analytic solution is available at present. For rota- 
tionally symmetric solutions, we write 



X|=/(r)e'"^ N = N{r), 



e^'^x^ 



9i(r), Ai ^ ^—[g2{r) - n], 



with various radial functions here satisfying the boundary conditions 



n/(0) 
^2(0) 



= yi(0) = iV'(0)-0, 



5(2(00) = h{oo) = 0, 



For these solutions, the angular momentum is given as 



'9^3 J 



1 



cfxe'^x' 



9'l{9l+92) + — 5-«/l^2 



fdr2n (^191 + 9192) 



0. 



(55) 



(56) 



(57) 



This null result is expected since there is no charged vortices in this model. 

In Fig. 4, some numerical results arc given for rotationally symmetric solutions 
with n — 1. Note that the profiles of \x\ and F12 are similar to those found in 
the ordinary abehan Higgs model. The Chern-Simons magnetic field /12 shows an 
interesting behavior. It is positive for small r, then changes sign and goes to as 
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r oo. This is expected because there is no charged vortex, i.e., = 0, as we 
noted above. 

The following comment on this model may be useful. One may view this model 
as the ordinary self-dual abelian Higgs model modified by the addition of the Chern- 
Simons term, and the effect is to make the neutral scalar nontrivial. But the main 
features of this model are essentially the same as those of the abelian-Higgs model: 
there is only broken phase, and no charged or spinning vortex exists. 

D. ei = 63, 62 = 64 

Here, two complex scalar fields x ^i-nd if) enter the theory in such a way that an 
additional global SU{2) symmetry may be present. The potential is 

U = ^(eal^P - u^) + ^1^1' feil^l' - k-n] , (58) 
2 V ei / 

where \l/ = ^ ^ ^ . The solution of this case will possess properties corresponding to 
a juxtaposition of the cases A and C; in particular, there is no charged or spinning 
vortex solution. As we have seen in the cases A and B, SU{2) global symmetry 
enables us to write w{z) = ^ = p"|^j , where P„, Qn are the same as before. Then 
the resulting self-duality equations can be read from Eqs. (p8| ) and (0), 

V^r/ - 2^e''(e'' - iV - 1) - e" + 1 = p, 

V^A^ - 6^(6" - - 1) = 0, (59) 

where 1] and p are defined by Eqs. ( p7|) and (^). In this case, however, a rigorous 
existence proof for p = case is not available and, strictly speaking, it is just a 
conjecture that the solution exists for each P„(z) and Qn{z) in Eq. (p^). 
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Without going further on this model, we just make one comment. As before we 
may consider this model as a linearized version of CP^ model to which a Chern- 
Simons gauge field is coupled. A neutral field A^ appears to make the system self- 
dual. From the self-duality equations (0), it is easy to see that for solutions of the 
form Eq. (H), 



(60) 



M ~ 0, 

A^ ~ when |g| ^ 1 
and they again approximate CP^ lumps as expected. 



E. ei = 63, 62 = —64 and u 



2 _ 







u 



Finally, we will consider the case with ei = 63 and 62 = —64, while setting 
= 0. This case is interesting because both gauge fields couple to scalar fields in a 



nontrivial manner, and there are both symmetric and asymmetric phases which are 
degenerate. The potential reads 



U 



Kixr-i^i 



2\2 



eidxr + lV^P 



ei 



1 2 



p2 



Ke2 



ei{\x\' + \m + —N 



ei 



(61) 
and 



and we have a symmetric phase if the vacuum corresponds to |x| = 14' 
N =(an arbitrary constant), and an asymmetric phase with the vacuum described 
by = I^P = v^/2ei and A^ = 0. 

To make the fields dimensionless, we make the replacements 



X,. 



N 



X 



X 



V 



2„,2 



-AT, 



V2el 
e2f 



3/2^/^' 



V 



V 
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Then after eliminating gauge fields from the self-duality equations we obtain the 
following set of equations: 

v^in \x? = e[(xr + m')ix\' + iv^r - 2)(ixr - m')N] + ixr - iv^r 

+47r5^5(x-x^), 

r=l 

VHn\^P\' = aixl' + l^n(xP + 2)(|xP - m')N] + IxP - \^\' (62) 

+47r^5(x-x',), 

r=l 

v^iv = -dxp - m'){x\' + i^r - 2) + ivdxr + i^r), 

where ^ = eff^/K^e^. Being nontrivially coupled and highly nonlinear, analytic 
solutions to these equations appear to be out of question. But with the ansatz 

1x1 = 1^1, iV = 0, (63) 

we can reduce the above three equations to the following single equation 

rii 

In \x\' - mxlMxl' - 1) = 4vr 5: 5(x - x,), 

r=l 

which is identical to that appearing in the self-dual Chern-Simons model with one 
scalar field [Q. (See also Eq. (0)-) It is unclear whether there exist solutions other 
than these. Especially, it will be interesting to see whether there is a solution with 
half integer vorticity for both and 0- 

Before investigating some of these issues, we first calculate quantum numbers of 
rotationally symmetric configurations, which are specified as 

X = /i(r)e^"^^ V^ = /2(r)e™^^ 

= e'-[g,{r) ^j, A, = e^-[g2{r) ^j, (64) 

N = N{r). 
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To ensure finite energy, various functions liere are subject to following boundary 
conditions: at spatial infinity, 



while, at r = 0, 



/i(oo) 
N{oo) 



/2(oo) = 1 or 0, 

a + (3 



92[oo) 



a — (3 



(65) 



^2/2(0) = iV'(0) = 0, 



ni + n2 



^?2(0) 



rii - n2 



(66) 



(67) 



[In Eq. (|65D, we must set a = j3 = Nq = {) for solutions in the broken phase.] Then, 
as usual, we find the flux values 



2 2 
while the angular momentum J is given by 



(6^ 



J 



2^62 J r r J 

27TV ' -2 



[A 

elea \2^ 
^ a + (3 

TTK, Ci 

ni + 722 
7CK, ei 



Q 
Q, 



(69) 
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where Q = Q^^ + Q^{=2Q^, from Eq. (^) denote the total electric charge. Note that 
J does not depend on couplings associated with the the Maxwell field; this is natural 
since Chern-Simons interaction is responsible for non-zero angular momentum. 

Finally we turn to self-duality equations in Eq. ( p^ ) and discuss the possibility of 
vortex solutions with half-integer vorticity (i.e. $a = 2ti{^+^y^) and PMa = 2vr(|+ 
^Y^))- Since Eq. (^) is very comphcated we made some numerical studies assuming 
the rotationally symmetric ansatz (|6^ ) with ni = 1, n2 = 0. Specifically we tried 



relaxation methods using the package COLSYS [15| as well as shooting methods 



both unidirectional and bidirectional [fH]. But we failed to find such solutions 
for both topological and nontopological cases. A plausible reason for our negative 
finding might be that the boundary condition A^(oo) may well be inconsistent with 
one of the Gauss laws which requires that = Q^. As we shoot the fields near the 
origin, we found that is either logarithmically divergent or becomes singular at a 
finite value of r. [Note that with the ansatz (|55D, = is automatically satisfied.] 
Of course, as we have already seen in one instance (at the end of the subsection A), 
that we have appropriate self-duality equations does not imply the actual existence 
of a non-trivial solution for which the given energy bound is saturated. In spite of 
what we have observed, however, there remains some possibility that such solutions 
do exist and our numerical analysis is incomplete; for example, solutions may be 
found with initial parameters which we did not sweep. Therefore, it is a bit hasty to 
conclude that there is no vortex solution with half-integer vorticities in this model. 

IV. SUMMARY 



In this paper we have studied various interesting special cases in the U{1) x U (1) 
self-dual Maxwell-Chern-Simons model. They include abelian Higgs system with two 
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complex scalar fields (the case A) and its Cliern-Simons counterpart (the case B). 
The latter can be regarded as a linearized version of the CP^ model with the Chern- 
Simons term (while maintaining self-duality) and as such it interpolates two familiar 
field theory systems appearing in the discussion of high- Tc superconductivity 
We proved the existence of multisoliton solutions for the topological soliton case, 
along the lines taken in Refs. IQ and ||^. In the symmetric phase we discussed 
the properties of nontopological soliton solutions by utilizing some known exact 
solutions of the corresponding nonrelativistic model. A more general field theory 
model interpolating the self-dual systems of the cases A and B is given in the 
appendix. 

With 63 = 64 = 0, the resulting self-dual model consists of a complex scalar, 
a neutral scalar, the Maxwell field, and the Chern-Simons field. It has no charged 
vortex because of the Gauss law and has more or less the same properties as those of 
the ordinary abelian Higgs model. The role of the Chern-Simons field is to modify 
field configurations such that neutral field may have nontrivial values. Extension 
of this case to a system with two complex scalar fields was also discussed. 

Finally we considered a model where all scalar fields and gauge fields couple 
to one another nontrivially. We calculated the quantum numbers for rotationally 
symmetric configurations, and showed that the ansatz |x| = leads back to the 
self-duality equations found in the 'minimal' self-dual Chern-Simons system [Q . We 
performed numerical analysis to find vortex solutions with half-integer vorticities 
and the result was negative. However, there is still a room for further examination. 

As we mentioned in the introduction, this generalized abelian Higgs model with 
independent Chern-Simons gauge term may have some relevance in describing real 
physical systems. Then a variety of vortex structures we discussed in this paper 



25 



are expected to have some useful physical applications. These issues are under 
investigation. 
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APPENDIX 

Here, we briefly describe the self-dual topologically massive U{1) gauge theory 
with two complex scalar fields and a neutral one, extending the model of Ref. such 
that the model may have an additional global SU (2) symmetry. The Lagrangian is 
given by 

£ = + \e^''^F,^A, + \D,^\' + ^{d,Nf - U{^, N), (70) 

where 

D^^ = {d^ - leA,,) 

and 

U{m, N) = -(e|^|2 + kAT - ev'^f + e^A^^l^j^ (71) 

with 

|^|'=|0l|'+|02p. 
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It has global SU (2) symmetry in addition to local U (1) symmetry. Gauss law reads 



diFio + kFi2 - eJ° = 0, 



(72) 



where 



Integrating over space, we obtain the usual relation 

where Q — JcfxJo. It is not difficult to show that the energy functional 
E = jd^x 



(73) 



(74) 



satisfies the following inequality 



E > lev^^l, 



(75) 



where the equality holds if and only if the following self-duality equations are satisfied 
(static fields assumed): 

{D^±iD2)^ = 0, 

Fio±diN = 0, 

(Do^fieN)^ = 0, (76) 

Fi2±{kN + e\^\^ -ev^) = 0. 

For a finite energy, we must here demand that 



r — > oo : 1^1 



v'^, N — > asymmetric phase. 
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or 

r — >■ oo : 1^1 — > 0, — ^ symmetric phase. 

K 

In the asymmetric phase, the vacuum manifold is 5''^. Therefore, as in the case A 
of Sec. Ill, the phase of the scalar field ^ fibers the vacuum manifold as a U{1) 
bundle over CP^ ~ S*^, and so a given sector of the theory is specified by a choice 
of a point on S"^ and the winding number of the scalar field around the fiber over 
that point. 

If only one complex scalar field were present, the model would interpolate the 
Ginzburg-Landau model and the pure Chern-Simons model (as noted in Ref. p). 
Hence, in the case of the above model (with SU{2) global symmetry), it should be 
clear that the case A and the case B of Sec. Ill can be obtained in the appropriate 
limiting cases. In other words, the solutions in this model will interpolates three 
kinds of solitons solutions, i.e., Ginzburg-Landau vortices, Chern-Simons vortices 
and CP^ lumps. 
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Figure captions 

Figure 1. Plot of solutions for the case A. (a) Rotationally symmetric solutions 
with n — 1, for p{= 64/62) = 1 (dotted lines) and for p = 2 (solid hne). We 
choose |x(oo)p = \xo\'^ = M^/2e2. Lines approaching zero asymptotically are 
magnetic fields F12. The solid line approaching 1 is that approaching l/\/2 
is IV'I, and dashed line approaching 1 is |x|(= for p — 1. (b) The functions 
6^/2 (approaching 1) and F12 (approaching 0) for rii — 1 and 712 = in the 
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62 = 64 = e case. The dotted, solid, dashed, and dot-dashed hnes correspond 
to |g| = 0.2, 1, 5, and 10, respectively. 

Figure 2. Plot of solutions in an asymmetric phase for the case B. (a) Rotationally 
symmetric solutions with n = 1 for p{= 63/61) = 1 (dotted lines) and for 
p = 2 (solid line). We choose |x(c)o)p = |xoP = v'^/2ei. Lines approaching 
zero asymptotically are magnetic fields /12. The solid line approaching 1 is 
that approaching l/\/2 is and the dashed line approaching 1 is |x|(= IV^I) 
for p = 1. (b) 6''/^ (approaching 1) and /12 (approaching 0) for solutions having 
the form in Eq. (|55D. Dotted, solid, dashed, and dot-dashed lines correspond 
to g = 0.2, 1, 5, and 10, respectively. 

Figure 3. Plot of solutions in the symmetric phase for the case B. (a) Solid lines are 
for n = and a = 2.1 with p = 63/61 = 2. Dashed lines are for n = 1 and 
a = 3 with p = 2. Among the same types of lines having values less than 1, 
upper lines correspond to |x| and the other to The magnetic field /12 is 
ring-shaped in both cases, (b) The functions 6^/^ (starting at ^^) and /12 for 
solutions of the form in Eq. (|39D. The dashed line is for |g| = 5, and solid line 
for |g| = 10. The dotted line is the plot of the function v^T2|g|/(r^ + |gp) for 
|g| = 10, a solution of nonrelativistic self-duality equations. 

Figure 4. Plot of rotationally symmetric solutions with unit vorticity for the case 
C. Dashed and solid lines are for ^ = 1 and 2, respectively. |x| are represented 
by lines approaching 1, and by negative definite lines. Lines starting at 1 
represent F12; lines having both positive and negative values represent /12. 
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